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Abstract. Experimental results obtained at the Relativistic Heavy Ion Collider 
(RHIC) have been interpreted in terms of a strongly interacting quark gluon 
plasma. The strongly interacting plasma is characterized by "perfect fluidity" , 
i.e. a ratio of shear viscosity to entropy density that saturates a proposed 
lower bound. In this contribution we explore the possibility that a similar 
phenomenon takes place in a strongly coupled non-relativistic Fermi liquid in 
^ ' which the scattering length between the Fermions is infinitely large. 

a; 

Keywords: hydrodynamics, viscosity, collective flow 
PACS: 25.75.Ld,66.20.+d 

> 

o 

■ 1. Introduction 



Experiments at the Relativistic Heavy Ion Collider (RHIC) indicate that a new 
C^, • state of matter is produced in high energy heavy ion collisions [ [T] . Much effort is 

, currently devoted to characterizing the properties of this state, and to determining 

the nature of its low energy excitations. 

Asymptotic freedom implies that the equation of state of a quark gluon plasma 
at T ^> Aqcd is that of a free gas of quarks and gluons. Numerical results from 
lattice QCD calculations show that at T ~ 2T C , which is relevant to the early stages 
of heavy ion collisions at RHIC, the pressure and energy density reach about 85% 
of the free gas limit. The 15% reduction is consistent with the magnitude of the 
first order perturbative correction. Higher order terms in the perturbative expan- 
sion converge very slowly, but this problem can be overcome using resummation 
techniques [ 2] . In this framework the degrees of freedom are dressed quasi-quarks 
and quasi-gluons, and these quasi-particles are weakly interacting. 

Transport properties of the matter created at RHIC indicate that this may 
not be correct. Experiments at RHIC suggest that the viscosity of the plasma is 
very small, and that the opacity for high energy jets is very large. If the plasma 
is composed of weakly interacting quasi-particles then the shear viscosity can be 
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estimated using kinetic theory. The result is 

V - ^npl, (1) 

where n is the density, p is the average momentum, and I is the mean free path. In 
a relativistic system the number of particles is not conserved, and it is more natural 
to express the viscosity in units of the entropy density rather than the density. Since 
the entropy per particle (in units of kg) is of order one, this does not qualitatively 
change the numerical coefficient in equ. ((T|). A good quasi-particle is characterized 
by a small ratio of the width over the excitation energy. This implies that the 
product of the mean free path times the typical momentum is large, and rj/s is big. 
This is confirmed by a weak coupling calculation in perturbative QCD. Arnold et 
al. find [[3] 

s 5 4 log(2.42 ff -i)' 1 1 

where g is the strong coupling constant. For a s < 1/3 we get rj/s > 1.75. This 
result should be contrasted with the values extracted at RHIC, which are in the 
range rj/s < 0.5 [0J M- 

From a theoretical point of view the RHIC results raise the question of how small 
the viscosity can get. Clearly, rj/s decreases as the interaction becomes stronger but 
there are good reasons to believe that the shear viscosity always remains finite. In 
particular, it seems reasonable to assume that the product of the mean free path and 
the typical momentum cannot be smaller than h [[6]. An interesting perspective on 
this issue is provided by a strong coupling calculation performed in the large N c limit 
of J\f = 4 SUSY Yang Mills theory. The calculation is based on the duality between 
strongly coupled gauge theory and weakly coupled string theory on AdSs x S§ 
discovered by Maldacena [ |7j . The thermodynamics of the SUSY field theory was 
studied by Gubser et al. [ [8] . They find that the entropy density in the strong 
coupling limit is 3/4 of the free field theory result. This implies that the equation 
of state is not drastically affected by the value of the coupling. The calculation was 
extended to transport properties by Policastro et al. These authors find that 
the shear viscosity to entropy density ratio of the strongly coupled gauge theory 
is rj/s = 7i/(47r). This number is quite consistent with the values extracted from 
RHIC data. Kovtun et al. studied the behavior of rj/s in other strongly coupled 
field theories with gravity duals and conjectured that the value K/ (4-7r) is a universal 
lower bound for rj/s [[TU1 ITTj . 

2. Cold atomic gases 

In order to understand the relevance of these results to the RHIC data it is useful to 
study the transport properties of other strongly coupled fluids that are experimen- 
tally accessible. Over the last ten years there has been truly remarkable progress 
in the study of cold, dilute gases of fermionic atoms in which the scattering length 
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a of the atoms can be controlled experimentally. These systems can be realized in 
the laboratory using Feshbach resonances, see [ Q2] for a review. A small negative 
scattering length corresponds to a weak attractive interaction between the atoms. 
This case is known as the BCS (Bardeen-Cooper-Schricffcr) limit. As the strength 
of the interaction increases the scattering length becomes larger. The scattering 
length diverges at the point where a bound state is formed. This is called the 
unitarity limit, because the scattering cross section saturates the s-wave unitarity 
bound a = Air/k 2 . On the other side of the resonance the scattering length is 
positive. In the BEC (Bose-Einstein condensation) limit the interaction is strongly 
attractive and the fermions form deeply bound molecules. 

The unitarity limit is of particular interest. In this limit the atoms form a 
strongly coupled quantum liquid which exhibits universal behavior. At zero tem- 
perature the atomic liquid is characterized by the mass of the atoms m, the density 
n, the scattering length a, and the effective range r. A dilute gas at unitarity corre- 
sponds to the limit in which a 3 n — > oo and r 3 n — ► 0. This implies that the depen- 
dence of a physical observable on n and m is determined by simple dimensional 
analysis, but the overall numerical constant is a complicated, non-perturbative 
quantity. At finite temperature T properties of the fluid are universal functions 
of the dimensionless variable T/Tp, where Tp ~ n 2 / 3 /m. 

In cold atomic gases we can reliably compute rj/s in the BCS limit [[13]. The 
ratio is temperature dependent and has a minimum at T ^ Tp. The shear viscosity 
is proportional to 1/a 2 , and rj/s is very large in the weak coupling limit. As in 
the case of QCD there are no controlled calculations in the strong coupling limit 
a — > oo. It is possible, however, to reliably extract rj/s from experimental data on 
the damping of collective oscillations [ Q31 HH [TB] . 

3. Collective Oscillations 

In the strong coupling limit we can assume that collective modes are approximately 
described by ideal fluid dynamics. The equation of continuity and of momentum 
conservation are given by 

0, (3) 

—VP - nW, (4) 

where n is the number density, m is the mass of the atoms, v is the fluid velocity, P 
is the pressure and V is the external potential. The trapping potential is harmonic 



- + v-M = 

dv {_ -a _ 

mn— + mn yy ■ V j v = 



(5) 
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Universality implies that the equation of state is given by P = rt 7+1 /(T/Tp) with 
7 = 2/3. The compressibility at constant entropy is 

dP\ P 

The equilibrium distribution hq can be determined from the hydrostatic equation 
VP = -n W. At T = 

no(0=no(0)(l-£-j|) , % = ^ (7) 

where /x is the chemical potential. In the unitarity limit the chemical potential is 
related to the Fermi energy as /i = £Ef, where £ is a universal parameter. Green 
function Monte Carlo calculations give £ ~ 0.44 [Hi]. We consider small oscillations 
n = tiq + <Jn. From the linearized continuity and Euler equation we get [ 118) 

m o^ = -7 (v • u) ( w) - v (tr • w) , (8) 

where we have dropped terms of the form VjVji? that involve higher derivatives of 
the velocity. This equation has simple scaling solutions of the form Vi = aiXi exp(ia;i) 
(no sum over i). Inserting this ansatz into equ. (jHJ we get an equation that deter- 
mines the eigenfrencies u). The experiments are performed using a trapping potential 
with axial symmetry, u>i = u>2 — u>o, = Xujq. In this case we find one solution 
with lj 2 = 2u} 2 and two solutions with [ [18l [HI HQ] 



- 2 =- 2 |7 + l + ^A 2 ±y^I±^A4 + (7 2 -37-2)A2 + (7 + l) 2 |. (9) 

In the unitarity limit (7 = 2/3) and for a very asymmetric trap, A — > 0, the 
eigenfrequencies are lu 2 = 2uj 2 and ui 2 = (10/3)o;q. The mode lj 2 — (10/3)wq is a 
radial breathing mode with a = (a,a,0) and the mode 10 2 = 2wq corresponds to a 
radial dipolc a = (a, —a, 0). 

The frequency of the radial breathing mode agrees very well with experimental 
results [ [3T] . In a hydrodynamic description the damping of collective modes is due 
to viscous effects. The dissipated energy is given by 

E=~J d 3 xTj(x) ^diVj + djVi - ^SijdkVkJ - J d?x(,(x) (dm) 2 , (10) 

where r\ is the shear viscosity and £ is the bulk viscosity. In the unitarity limit the 
system is scale invariant and the bulk viscosity in the normal phase vanishes. For 
the radial scaling solutions we get 

E = -| (4 +a 2 y - a x a y ) j d 3 xr)(x), (11) 
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Fig. 1. Viscosity to entropy density ratio of a cold atomic gas in the unitarity 
limit, from [[H]. This plot is based on the damping data published in and the 
thermodynamic data in [[2H H5]. The dashed line shows the conjectured viscosity 
bound rj/s = 1/(Att). 



where E is a time average. The damping rate is given by the ratio of the energy 
dissipated to the total energy of the collective mode. The kinetic energy is 



m f , „ j0 mN 
~2 



E km = - I d 3 xn(x)v 2 = — {al + a 2 ) (x 2 ). (12) 



At T — we find (x 2 ) — i?j_/8, where R± is the transverse size of the cloud. At 
non-zero temperature we can use the Virial theorem [[55] to relate (x 2 ) to the total 
energy of the equilibrium state, (x 2 ) / {x 2 )t=o = E/Et=o- The damping rate is [ 

nans] 

I /'.' 2 a 2 + a 2 - a x a y fd 3 xrj(x) 



2E 3 al + a 2 mN(x 2 



(13) 



Note that the second factor on the RHS is 1/2 for the radial breathing mode and 
3/2 for the radial dipole mode. If this dependence could be demonstrated experi- 
mentally, it would confirm that the damping is indeed dominated by shear stress. 

We shall assume that the shear viscosity is proportional to the entropy density, 
rj(x) = as(x). We note that since the flow profile has a simple scaling form the 
damping rate is proportional to the volume integral of the shear viscosity. If 77 ~ s 
then the damping rate is proportional to the total entropy. The kinetic energy, 
on the other hand, scales with the number of particles. We can now relate the 
dimensionless damping rate T/luj_ = l/(ru±) of the radial breathing mode to the 
shear viscosity to entropy density ratio. We find 
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Fig. [T] shows rj/s extracted from the experimental results of the Duke group [I23j. 
The entropy per particle was also taken from experiment [I25j. Similar results are 
obtained if the entropy is extracted from quantum Monte Carlo data. The critical 
temperature for the superfluid/normal transition is T c /Tp — 0.29. We observe that 
rj/s in this regime is roughly 1/2. This value is compatible with the conjectured 
viscosity bound and comparable to the values that have been extracted at RHIC. 



4. Elliptic Flow 

Collective modes are very useful because it is possible to track many compression 
and expansion cycles and even small damping coefficients can be measured. In 
heavy ion collisions we do not have this luxury and we have to rely on collective flow 
measurements to extract transport coefficients. In the following we shall estimate 
the effect of a non-zero shear viscosity on the elliptic flow of a cold atomic liquid. We 
consider the expansion of the atomic cloud after the trapping potential is removed. 
The expansion is described by a simple scaling ansatz 

bAt)ri 

n(ri,t) = n (ri/bi(t)) (i = l,...,3), Vi(f,t) = - 1 . (15) 

It is easy to check that this ansatz satisfies the continuity equation. The Euler 
equation gives [[55] 

k - (16 » 

We are interested in an axially symmetric trap with lu x = u) y = uj±_ and lu z = Xuj±. 
In the limit b± ^> b z we get 

L-^, S.-£ (17) 

The solution of these equations for a very asymmetric trap (A = 0.045) is shown 
in Fig. O We observe the usual elliptic flow phenomenon: the transverse pressure 
exceeds the longitudinal pressure, there is more acceleration in the transverse direc- 
tion, and as a consequence the transverse expansion is faster than the longitudinal 
one. The right panel shows that during the expansion internal energy is converted 
to kinetic energy of the flow. We can also compute the amount of energy dissipated 
due to viscous effects. We find 




d 3 x r)(x). (18) 

An estimate of this quantity is also shown in the right panel of Fig. [2] We show two 
curves, corresponding to rj/s equal to one and three times the conjectured viscosity 
bound, respectively. We have taken the entropy s at T = T c . We observe that most 
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Fig. 2. The left panel shows the evolution of the scale factors b± and b z as a function 
of the dimensionless variable Lu±t. The right panel shows the time evolution of the 
internal and kinetic energies. We also show an estimated of the energy dissipated 
by viscous effects. The two curves correspond to rj/s equal to one and three times 
the conjectured viscosity bound (using s at T = T c ). 



of the energy is dissipated early during the expansion. The total energy dissipated 
amounts to a few percent of the total energy available. This should make the effect 
observable, although it seems unlikely that a measurement of rj/s based on elliptic 
flow can be as accurate as the one based on collective modes. A measurement of 
elliptic flow was reported in [ 127] , but there is no systematic study of the temperature 
dependence of the effect. 
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